Contrasting patterns of population structure and gene flow facilitate exploration of connectivity in two widely distributed temperate octocorals

This is the final version of the article. Available from Springer Nature via the DOI in this record.Connectivity is an important component of metapopulation dynamics in marine systems and can influence population persistence, migration rates and conservation decisions associated with Marine Protected Areas (MPAs). In this study, we compared the genetic diversity, gene flow and population structure of two octocoral species, Eunicella verrucosa and Alcyonium digitatum, in the northeast Atlantic (ranging from the northwest of Ireland and the southern North Sea, to southern Portugal), using two panels of thirteen and eight microsatellite loci, respectively. Our results identified regional genetic structure in E. verrucosa partitioned between populations from southern Portugal, northwest Ireland, and Britain/France; subsequent hierarchical analysis of population structure also indicated reduced gene flow between southwest Britain and northwest France. However, over a similar geographical area, A. digitatum showed little evidence of population structure, suggesting high gene flow and/or a large effective population size; indeed, the only significant genetic differentiation detected in A. digitatum occurred between North Sea samples and those from the English Channel/northeast Atlantic. In both species the vast majority of gene flow originated from sample sites within regions, with populations in southwest Britain being the predominant source of contemporary exogenous genetic variants for the populations studied. Unsurprisingly, historical patterns of gene flow appeared more complex, though again southwest Britain appeared an important source of genetic variation for both species. Our findings have major conservation implications, particularly for E. verrucosa, a protected species in UK waters and listed by the IUCN as ‘Vulnerable’, and for the designation and management of European MPAs.We thank Natural England (project No. RP0286, contract No. SAE 03-02-146), the NERC (grant No. NE/L002434/1) and the University of Exeter for funding this research. Additional funding for sample collection, travel and microsatellite development was provided by the EU Framework 7 ASSEMBLE programme, agreement no. 227799, and NERC grant No. NBAF-362

Electron Transfer Function versus Oxygen Delivery: A Comparative Study for Several Hexacoordinated Globins Across the Animal Kingdom

Caenorhabditis elegans globin GLB-26 (expressed from gene T22C1.2) has been studied in comparison with human neuroglobin (Ngb) and cytoglobin (Cygb) for its electron transfer properties. GLB-26 exhibits no reversible binding for O2 and a relatively low CO affinity compared to myoglobin-like globins. These differences arise from its mechanism of gaseous ligand binding since the heme iron of GLB-26 is strongly hexacoordinated in the absence of external ligands; the replacement of this internal ligand, probably the E7 distal histidine, is required before binding of CO or O2 as for Ngb and Cygb. Interestingly the ferrous bis-histidyl GLB-26 and Ngb, another strongly hexacoordinated globin, can transfer an electron to cytochrome c (Cyt-c) at a high bimolecular rate, comparable to those of inter-protein electron transfer in mitochondria. In addition, GLB-26 displays an unexpectedly rapid oxidation of the ferrous His-Fe-His complex without O2 actually binding to the iron atom, since the heme is oxidized by O2 faster than the time for distal histidine dissociation. These efficient mechanisms for electron transfer could indicate a family of hexacoordinated globin which are functionally different from that of pentacoordinated globins

Population genomics applications for conservation: the case of the tropical dry forest dweller Peromyscus melanophrys

Recent advances in genomic sequencing have opened new horizons in the study of population genetics and evolution in non-model organisms. However, very few population genomic studies have been performed on wild mammals to understand how the landscape affects the genetic structure of populations, useful information for the conservation of biodiversity. Here, we applied a genomic approach to evaluate the relationship between habitat features and genetic patterns at spatial and temporal scales in an endangered ecosystem, the Tropical Dry Forest (TDF). We studied populations of the Plateau deer mouse Peromyscus melanophrys to analyse its genomic diversity and structure in a TDF protected area in the Huautla Mountain Range (HMR), Mexico based on 8,209 SNPs obtained through Genotyping-by-Sequencing. At a spatial scale, we found a significant signature of isolation-by-distance, few significant differences in genetic diversity indices among study sites, and no significant differences between habitats with different levels of human perturbation. At a temporal scale, while genetic diversity levels fluctuated significantly over time, neither seasonality nor disturbance levels had a significant effect. Also, outlier analysis revealed loci potentially under selection. Our results suggest that the population genetics of P. melanophrys may be little impacted by anthropogenic disturbances, or by natural spatial and temporal habitat heterogeneity in our study area. The genome-wide approach adopted here provides data of value for conservation planning, and a baseline to be used as a reference for future studies on the effects of habitat fragmentation and seasonality in the HMR and in TDF

Local selection in the presence of high levels of gene flow: Evidence of heterogeneous insecticide selection pressure across Ugandan Culex quinquefasciatus populations

Background: Culex quinquefasciatus collected in Uganda, where no vector control interventions directly targeting this species have been conducted, was used as a model to determine if it is possible to detect heterogeneities in selection pressure driven by insecticide application targeting other insect species. Methodology/Principal findings: Population genetic structure was assessed through microsatellite analysis, and the impact of insecticide pressure by genotyping two target-site mutations, Vgsc-1014F of the voltage-gated sodium channel target of pyrethroid and DDT insecticides, and Ace1-119S of the acetylcholinesterase gene, target of carbamate and organophosphate insecticides. No significant differences in genetic diversity were observed among populations by microsatellite markers with HE ranging from 0.597 to 0.612 and low, but significant, genetic differentiation among populations (FST = 0.019, P = 0.001). By contrast, the insecticide-resistance markers display heterogeneous allelic distributions with significant differences detected between Central Ugandan (urban) populations relative to Eastern and Southwestern (rural) populations. In the central region, a frequency of 62% for Vgsc-1014F, and 32% for the Ace1-119S resistant allele were observed. Conversely, in both Eastern and Southwestern regions the Vgsc-1014F alleles were close to fixation, whilst Ace1-119S allele frequency was 12% (although frequencies may be underestimated due to copy number variation at both loci). Conclusions/Significance: Taken together, the microsatellite and both insecticide resistance target-site markers provide evidence that in the face of intense gene flow among populations, disjunction in resistance frequencies arise due to intense local selection pressures despite an absence of insecticidal control interventions targeting Culex

Genetic Structure, Linkage Disequilibrium and Signature of Selection in Sorghum: Lessons from Physically Anchored DArT Markers

Population structure, extent of linkage disequilibrium (LD) as well as signatures of selection were investigated in sorghum using a core sample representative of worldwide diversity. A total of 177 accessions were genotyped with 1122 informative physically anchored DArT markers. The properties of DArTs to describe sorghum genetic structure were compared to those of SSRs and of previously published RFLP markers. Model-based (STRUCTURE software) and Neighbor-Joining diversity analyses led to the identification of 6 groups and confirmed previous evolutionary hypotheses. Results were globally consistent between the different marker systems. However, DArTs appeared more robust in terms of data resolution and bayesian group assignment. Whole genome linkage disequilibrium as measured by mean r2 decreased from 0.18 (between 0 to 10 kb) to 0.03 (between 100 kb to 1 Mb), stabilizing at 0.03 after 1 Mb. Effects on LD estimations of sample size and genetic structure were tested using i. random sampling, ii. the Maximum Length SubTree algorithm (MLST), and iii. structure groups. Optimizing population composition by the MLST reduced the biases in small samples and seemed to be an efficient way of selecting samples to make the best use of LD as a genome mapping approach in structured populations. These results also suggested that more than 100,000 markers may be required to perform genome-wide association studies in collections covering worldwide sorghum diversity. Analysis of DArT markers differentiation between the identified genetic groups pointed out outlier loci potentially linked to genes controlling traits of interest, including disease resistance genes for which evidence of selection had already been reported. In addition, evidence of selection near a homologous locus of FAR1 concurred with sorghum phenotypic diversity for sensitivity to photoperiod

Genetic Population Structure in the Antarctic Benthos: Insights from the Widespread Amphipod, Orchomenella franklini

Currently there is very limited understanding of genetic population structure in the Antarctic benthos. We conducted one of the first studies of microsatellite variation in an Antarctic benthic invertebrate, using the ubiquitous amphipod Orchomenella franklini (Walker, 1903). Seven microsatellite loci were used to assess genetic structure on three spatial scales: sites (100 s of metres), locations (1–10 kilometres) and regions (1000 s of kilometres) sampled in East Antarctica at Casey and Davis stations. Considerable genetic diversity was revealed, which varied between the two regions and also between polluted and unpolluted sites. Genetic differentiation among all populations was highly significant (FST = 0.086, RST = 0.139, p<0.001) consistent with the brooding mode of development in O. franklini. Hierarchical AMOVA revealed that the majority of the genetic subdivision occurred across the largest geographical scale, with Nem≈1 suggesting insufficient gene flow to prevent independent evolution of the two regions, i.e., Casey and Davis are effectively isolated. Isolation by distance was detected at smaller scales and indicates that gene flow in O. franklini occurs primarily through stepping-stone dispersal. Three of the microsatellite loci showed signs of selection, providing evidence that localised adaptation may occur within the Antarctic benthos. These results provide insights into processes of speciation in Antarctic brooders, and will help inform the design of spatial management initiatives recently endorsed for the Antarctic benthos

High Risks of Losing Genetic Diversity in an Endemic Mauritian Gecko: Implications for Conservation

Genetic structure can be a consequence of recent population fragmentation and isolation, or a remnant of historical localised adaptation. This poses a challenge for conservationists since misinterpreting patterns of genetic structure may lead to inappropriate management. Of 17 species of reptile originally found in Mauritius, only five survive on the main island. One of these, Phelsuma guimbeaui (lowland forest day gecko), is now restricted to 30 small isolated subpopulations following severe forest fragmentation and isolation due to human colonisation. We used 20 microsatellites in ten subpopulations and two mitochondrial DNA (mtDNA) markers in 13 subpopulations to: (i) assess genetic diversity, population structure and genetic differentiation of subpopulations; (ii) estimate effective population sizes and migration rates of subpopulations; and (iii) examine the phylogenetic relationships of haplotypes found in different subpopulations. Microsatellite data revealed significant population structure with high levels of genetic diversity and isolation by distance, substantial genetic differentiation and no migration between most subpopulations. MtDNA, however, showed no evidence of population structure, indicating that there was once a genetically panmictic population. Effective population sizes of ten subpopulations, based on microsatellite markers, were small, ranging from 44 to 167. Simulations suggested that the chance of survival and allelic diversity of some subpopulations will decrease dramatically over the next 50 years if no migration occurs. Our DNA-based evidence reveals an urgent need for a management plan for the conservation of P. guimbeaui. We identified 18 threatened and 12 viable subpopulations and discuss a range of management options that include translocation of threatened subpopulations to retain maximum allelic diversity, and habitat restoration and assisted migration to decrease genetic erosion and inbreeding for the viable subpopulations

Sex change and effective population size : implications for population genetic studies in marine fish

Large variance in reproductive success is the primary factor that reduces effective population size (Ne) in natural populations. In sequentially hermaphroditic (‘sex-changing’) fish, the sex ratio is typically skewed and biased toward the ‘first’ sex, while reproductive success increases considerably after sex change. Therefore, sex-changing fish populations are theoretically expected to have lower Ne than gonochorists (separate sexes), assuming all other parameters are essentially equal. In this study, we estimate Ne from genetic data collected from two ecologically similar species living along the eastern coast of South Africa: one gonochoristic, the ‘santer’ sea bream Cheimerius nufar, and one protogynous (female-first) sex-changer, the ‘slinger’ sea bream Chrysoblephus puniceus. For both species, no evidence of geneticstructuring, nor significant variation in genetic diversity, were found in the study area. Estimates of contemporary Ne were significantly lower in the protogynous species, but the same pattern was not apparent over historical timescales. Overall, our results show that sequential hermaphroditism may affect Ne differently over varying time frames, and that demographic signatures inferred from genetic markers with different inheritance modes also need to be interpreted cautiously, in relation to sex-changing lifehistories